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Abstract 

An S-matrix is proposed for the two dimensional 0(3) cr-model with a dynamical 
^-term (axion model). Exploiting an Abelian T-duality transformation connecting 
the axion model to an integrable SU(2)xU(l) symmetric principal a-model, strong 
evidence is presented for the correctness of the proposed S'-matrix by comparing the 
perturbatively calculated free energies with the ones based on the Thermodynamical 
Bethe Ansatz. This T-duality transformation also leads to a new Lax-pair for both 
models. The quantum non-integrability of the 0(3) a-model with a constant 9- 
term, in contradistinction to the axion model, is illustrated by calculating the 2^3 
particle production amplitude to lowest order in 9. 



In this paper we shall study the following two dimensional cr-model described by 
the LagrangianQ: 



where n"n" = 1. In fact (|I|) is an 0(3) non-linear cr-model, coupled to a scalar field, 6, 
(whose normalization has been chosen for later convenience) through the Hopf term. 
This latter is proportional to the topological current of the 0(3) model and with 
the normalization chosen in (|l]) its space-time integral (after a Wick rotation) yields 
the topological charge, which can take integer values only. This implies that the 
variable 6 is actually an angle, taking its values between and 2tt. This observation 
will play an important role in all our considerations. 

We shall refer to (|l]) as the 'axion model' since it can be thought of as the 
0(3) nonlinear a-model with a dynamical ^-term which can be regarded as a two- 
dimensional analogue of its phenomenologically important four-dimensional coun- 
terpart The following (heuristic) consideration might be useful to gain some 
insight into the physics of the axion model. Let us integrate out the 0(3) fields, 
n", in some generating functional of the theory (|l]). This way one would obtain a 
non-vanishing effective potential for the 6 field. Because 6 is 27r-periodic the effec- 
tive potential must be also periodic. The effective theory of the 6 field is therefore 
expected to be similar to the Sine-Gordon model, with a periodic potential and 
corresponding topological current = e^^d^O /2ti . Note that 6 being an angular 
variable makes it difficult to integrate it out in the functional integral in spite of C 
being only quadractic in 6. 

The axion model belongs to a family of (classically) integrable two-dimensional 
non-linear cr-models with an 0(3) symmetry discovered in Ref. 0. All of the models 
of Ref. [0] have target spaces with non-vanishing torsion (in addition to the metric 
tensor field) but the axion model is especially simple, as its torsion is constant. 

It has already been observed in Ref. that the axion model (|T|) can be mapped 
to a one parametric deformation of the SU(2) principal a-model by an Abelian 
T-duality transformation. This latter (torsionless) model has an SUl(2)®Ur(1) 
symmetry and recently there has been some revival of interest in it 0. Its 
Lagrangian can be written as: 

C^ = -^{LlL'^>' + gLlL'^}, (2) 

where 

L, = G-'d,G = T^Ll, (3) 

"'^We use the following conventions. For a vector v in two-dimensional Minkowski space v^v^ — 
Vj^v^ where v± = vq ±vi. The antisymmetric tensor is defined by e'^^ = 1, = with 
being the standard Pauli-matrices satisfying a'^a'' = 6°''' + ie'^^'^a'^. 
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and g is the parameter of deformation. The theory is known to be integrable 
classically p and there is little doubt that it is also quant um-integr able 0, 0. 
Its spectrum contains two massive doublets (kinks) whose scattering is described by 
the tensor product of an SU(2)xU(l) symmetric solution of the bootstrap S'-matrix 
equations: 

si^e) = s^'^\e) ® s^p\e) , (4) 

where S^''\9) denotes the Sine-Gordon (SG) S'-matrix. Depending on the value 
of the parameter, p, in addition to the kinks there are also some bound states 
(breathers) in the spectrum transforming as 3 + 1 under SU(2). 

In the following we shall show that the somewhat unexpected identification be- 
tween the axion model and the deformed principal cr-model through a T-duality 
transformation allows us to learn more about both of them. Assuming the validity 
of the duality transformation at the quantum level between the two theories implies 
the absence of particle production in the axion model (|l]), and also that its factorized 
scattering theory is given by the S'-matrix of Eq. (^. 

The proposed quantum integrability of the axion model might seem somewhat 
surprising, as it is generally believed that the 0(3) model with a constant ^-term 
is not quantum integrable, except for the special value 6^ = vr (despite the fact 
that the ^-term, being a total derivative, does not change the classical physics of the 
model). We now show that in the framework of the form-factor bootstrap approach 
the 6 term mediates particle production in the 0(3) cr-model, indeed. To lowest 
order in 6 the 2^3 particle production amplitude can be written as 

(p, 6; p\ y- f, b"\q, a; q', a')^e) = (2vr)2 1 9 6^'^ {p + p' + f - q - q') 

■ {p, b; p', b'; p", b"\T{0)\q, a; g', a')(o) + 0(9') , (5) 

where in the first line the amplitude is in the 0(3) model with a ^-term, while in the 
second line the matrix element of the topological charge density operator T is to be 
calculated in the original 0(3) cr-model (with = 0). In other words we simply apply 
perturbation theory in 9. Let us now consider the following simplified kinematical 
configuration: the incoming particles have momenta qi = Q and q[ = —Q, whereas 
the produced (outgoing) three particles have momenta pi = Q', p[ = and p" = —Q' 
respectively. Here Q' can easily be expressed in terms of Q and the kink mass M 
using energy conservation. For large Q, using the results of Ref. 0, we find 

{p,b;p\b';p'\b''\Tmq,a;q\a')^,^^nl^^^j^ . (6) 

Eq. (P) shows that already to first order in 9, the 2^3 particle production ampli- 
tude is different from zero. Thus at least for small values of 9, the introduction of 
this term destroys quantum integrability of the 0(3) a-model, indeed. 

To exhibit now the classical T-duality transformation between the two models 
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1^, we introduce the parametrization 

ra^ = sin "(9 sin (y9 , ra^ = sin -i? cos , = cos-d ^ 9 = — + (7) 

A 

in terms of which the Lagrangian (|1|) (after an integration by parts) becomes 

£ = i/^^^^M^ + sin^ §d^^d^^ + d^xd^x + 2^Jl + gcos^e^'^'d^xduv] ■ (8) 
2A 



We now perform an Abehan T-duahty transformation [10| with respect to the vari- 



able Xy which corresponds to the canonical transformation |11 



X = Tf^P" Px = r « ) (9) 

y ^ + 9 A 

where (and in the following) resp. pa denote the canonical momenta conjugate 
to X resp. to its 'dual' a. In terms of these new variables the dual Lagrangian turns 
out to be: 

£2 = ^{d^,dd^§+{l+gcoi^^d)d^,^d''^+{l+g) [d^ad^a + 2coi^dd^,ad^'^] } , (10) 
2A 

which is nothing but the Lagrangian when parametrizing the SU(2) valued field, 
G, by the Euler angles 

G = e'""^' e'^^' e^""' , (11) 
and taking into account the relations at the classical level between the couplings: 

A = A, g = g. (12) 

The observation, that the axionic model is the T dual of £s also explains why 9 
is an angular variable. Indeed it has been shown in Ref. [|12| that in case of the 
principal cr-model {g = 0) the Abelian T duality (||) maps its target space (S*^) into 
S'^ X S^. The arguments of Ref. [|l2l can be easily applied to the present case with 
g > —1, and it is clear that in Eq. (|1]) n"" parametrize the S"^ and 9 parametrizes the 
S\ 

In fact the equations of motion of both models (HD and (|]) are known to admit a 
Lax representation indicating their (classical) integrability f^, |^ . Indeed, introduc- 
ing the matrix valued current 

= ■^ne^^d'"9 - ^e^y&'n + \n^^Jl , (13) 

where n = in°'a"', the equations of motion of (|^) can be written as: 

a% = 0, d^I,~dJ,= [I„Q. (14) 
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The standard form (0) of the equations of motion allows for the introduction of a 
Lax pair 

U.-j^I., (15) 

satisfying the zero curvature equation 

d^U,-dMf.= [U^,U,], (16) 

for all values of the spectral parameter u. The current, is closely related to 
the matrix valued Noether current, A/"^ = —ir'^N'^, defined by 5C = d^e°'N'^ corre- 
sponding to the symmetry transformation = e^^'^e^n^: 

= W^ + e^.a^T, T=(^6-^)n. (17) 

The fact that, apart from a trivially conserved piece, can be identified with the 
Noether current of the manifest 0(3) symmetry of the Lagrangian explains only the 
first equation in Eq. (p!^. The trivially conserved part of is essential that the 
zero curvature equation be also satisfied. 

The equations of motion of the deformed principal model (0) can be written 
entirely in terms of the current as 



It is known that this system can be put to the Lax form |^ |^, i.e. there is a spectral 
parameter dependent current, = t^V^, satisfying the zero curvature equation 



(p^). This current can be written as: 

Vl'' = aiL'f, Vi = a^Ll, (19) 

where 

"^-~4-^^2±4^' «±--4-^a;2±4a; " 
We can now use the classical T-duality transformation M) to map the linear 



system of the axion model ([1^) to a new Lax pair for the deformed a-model (|]). It 



is given by Eq. (|l^), where the current, J^, has to be replaced by 

i, = d,G + g{Gr''G-')Ll - t^e^^ff" {Gt'G-') . (21) 



Eq. ( piD is obtained from (|T3D by the T-duality transformation @. is related to 
the Noether current TV^, corresponding to the manifest symmetry 6G = —ie"'T°'G of 
(H) and can be written analogously to J\f^: 

i,, = }Sf^ + e^,d''f, f = -i^GT'G-\ (22) 
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It is clear that the new Lax pair (|I5]) and the 'old' one, (0), cannot be related 
by a gauge transformation since they have different pole structures as functions of 
the spectral variable, uj. In the g —>■ —1 limit, the axion model reduces to the 
original 0(3) a-model (decoupled from the 6 field), and the Lax pair ([T5|) becomes 
equivalent to that of Ref . JTSf , where it has been pointed out that the corresponding 
J^'s are ultralocal currents. We note that the Lax pairs (|T5|) and (|19D correspond 
to (different) deformations of the usual Lax pairs of the principal chiral cr-model, 
linear in d^GG~^ respectively G^^d^G. 

Next we carry out a standard test on the proposed S'-matrix (^) for the axion 
model by comparing its (zero temperature) free energy obtained from the Thermo- 
dynamical Bethe Ansatz (TBA) and in weak coupling perturbation theory (PT) 
iM. For the deformed cr-model (0) this comparison has been done in Ref. [0 where 



complete consistency has been found between the results of PT and of the TBA. For 
the axion model it is sufficient to compute the free energy in PT as the results of 
the TBA can be literally taken over from Ref. 0. In the present case one obtains as 
a bonus, a further nontrivial check on the quantum equivalence between the axion 
and the deformed a-model, hence also on the validity of the T-duality transforma- 
tion at the quantum level. Up to now when quantum equivalence between dually 
related models has been tested, mostly /9-functions have been compared. The fact 
that the higher coefficients of the /3-functions are scheme dependent makes such a 
comparison more difficult and less conclusive. 

The equivalence of the /^-functions is certainly a necessary condition for the 
validity of quantum T-duality. At one loop order the /3-functions of the couplings, 
Px, 13 g and (3^^, [3g are simply obtained from each other by the classical relation 
(|l^). At two loops, however, it has been found in that using the background 
field method and dimensional regularization the following perturbative redefinition 
of the couplings 

A = A+^(l + (7), ~g=g + ^[l+g)\ (23) 

(i.e. a change of scheme) is induced by the T-duality transformation. Taking into 
account Eqs. (p3D the two loop /^-functions of the two models turn out to be equiva- 
lent. Alternatively, introducing a renormalization group (RG) invariant combination 
of the two couplings: 

p = 27r lim ' , p = 27r lim T ' , 24 

where t ccXnh, one finds p = p up to two loops . It is important to note that 
the RG invariant quantity ( p4|) can be consistently identified with the parameter p 
in the S matrix (^. Let us introduce an effective /5-function for X(t) by i3^ft{X,p) = 
/3a(A, r(A,p)), expressing g{t), in terms of the running coupling, X{t), and p as 
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g{t) = r(A(t) ,p). Using the perturbative result for T{X{t),p) p one finds 

/3.«(A,p) = A«(A,p) = + + ■■■■ (25) 

Tlius as far as coupling constant renormalization is concerned, the two models are 
equivalent, both are asymptotically free, and the actual value of p effects only the 
two loop coefficient. 

The classical free energy density is obtained by minimizing the Legendre trans- 
form of the Hamiltonian density coupled to some conserved currents 

n = no-hi4, H = J dxn = H -hiQi. (26) 

Since the axion field, 6, is actually an angle, its winding number (topological charge) 
can be non trivial. Therefore we present here the Legendre transformation of the 
modified Hamiltonian ( p6D for a rather general case. 

Let us consider a general sigma model with torsion 

Co = ^9ABd''X%X'' + h^ABe^'d^X^d^X'', (27) 

and the following Ansatz for a set of conserved currents 

j; = CAX)d,X^ + e^^B\{X)d^X^ , (28) 

sufficiently general to include topological currents. The Legendre transformation of 
(^) yields the Lagrangian of the modified model which can be written as 

t = Co + h'4 + hi'h^Cifi^^ . (29) 

In fact C can be obtained by gauging Cq i.e. by the substitution 

d^X^ d,X^ + h'5,,C^ (30) 

when the antisymmetric field hAs is invariant (without compensating gauge trans- 
formation) under the symmetry transformation generated by the conserved currents 



Below we also give a class of classical ground states (around which C is to be ex- 
panded) assuming that the metric, the antisymmetric tensor field and the quantities 
characterizing the currents are independent of a set of coordinates, 6'°, corresponding 
to the splitting X^ = {y'', ^"): 

gAB = 9AB{y), bAB = bAB{y), C\ = C\{y), B\ = B\{y). (31) 
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In this case the ground state is characterized by constant y'^-s and constant 6'"-s 
y'' = y^, 6'" = 6'q, where the y^-s stand for the extrema of 

H,s = -h^h^igABC^'C"^ + Bi^&p{r^r^) , (32) 



and 



^'o = {l-'r^BUyo)h\ (33) 



In ( |52| - P3D 7q,/3 denotes the restriction of qab to the submanifold coordinatized by 

The axion model has a 'manifest' (i.e. up to a total derivative) SU(2) xU6)(l) 
symmetry, where the U6i(l) subgroup is generated by the shift 9^9 + const. It 
is very important to note that although the SU(2) symmetry of the Lagrangian 
(|I]) corresponds to that of the S'-matrix, the 'manifest' U6i(l) symmetry cannot 
be identified with the corresponding one of the S'-matrix (^). Here the duality 
transformation provides the clue; the corresponding Ur(1) symmetry of the axion 
model is actually the image of the manifest Ur(1) symmetry of (0) under the duality 
transformation, thus it is generated by the topological current of the axion field. 

Corresponding to the SUl(2)(8)Upj,(1) symmetry of the deformed a-model there 
are two Noether charges, Ql resp. Qr, associated to the Ul(1) resp. Ur,(1) sub- 
groups. Introducing two chemical potentials coupled to the Ql resp. Qr, charges 



the Hamiltonian (26) takes the form: H = Hy. — /^lQl — hjiQ^ . Then one can 
distinguish between three different types of finite density ground states: LEFT with 
/iL > 0, and /iR = 0, RIGHT with /il = 0, and /ir > 0, and DIAG where /il, /^r > 0. 
As found in ||^ the RIGHT case is obtained from DIAG by letting /il = in the 
Gnal results. 

We compute below the corresponding ground state energies to one loop order 
in the axion model starting with the LEFT case first. With the Euler angle 
parametrization of G (|Tl]) the Ul(1) transformation, G t— > e*'*^^G of the deformed 
cr-model (0) acts as a simple shift, <f{x) ^ f{x) + k. The corresponding Noether 
charge, Ql, and its image under the T-duality transformation, Ql, are simply 



f ~ f djCY; dC 

Ql = J dxp^ , = J d-xp^ , where = ' -^"^ ^ ^ 



since the canonical transformation implementing the T-duality mapping (|^) effects 
only Pa, x' 1 ct' and p^, leaving the other fields, 'd, p^, p^, unchanged. 

Since in the LEFT case the bAs field in Eq. (|1|) is invariant, one can simply 
'gauge' the Lagrangian of the axion model in an external (/il) field (see Eq. (pO]). 
The classical ground state is found tobev9 = x = 0,'i9 = 7r/2. (The corresponding 
solution of the deformed cr-model is given hj ip = a = 0, = 7r/2.) 
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Expanding the (Euchdean) Lagrangian (after suitable rescalings, etc.) we obtain 



C 



2h 



+ -mMm^ + o(A), 
Ao ^ 



where 



M 



-d^ + Ahl 










(34) 



(35) 



and m = (t?, ip, x)- (Aq, Qo denote the bare coupling and parameter of the axion/dual 
model). In Eq. (^) we kept the e tensor explicitly, as it requires a careful definition 
in n = 2 — e dimensions which we use to regularize the momentum integrals. We 



adopt the definiton of [0 , where this antisymmetric tensor corresponds to an almost 
complex structure: e 



Sua- The one loop quantum corrections to 
the classical ground state (the first term in Eq. (|3^)) require the calculation of a 
functional determinant, leading to 



n 



(T-p 



(27r 



QoVl) 



(36) 



To evaluate (^) we apply the modified dimensional regularization of Ref. |^ as 
p2 = ^2vPu plays here a distinguished role and it is kept as a one dimensional variable. 
In fact for our purposes it is sufficient to calculate the difference J^{h) —J-'Y:{h), where 
^T,{h) is the corresponding determinant in the deformed cr-model (Eq. (3.12) in Ref. 
1^). Since both J-'{h) and TT.{h) are already the first quantum corrections to the 
classical expressions we may set g = g (and make no distinction between bare and 
renormalized g^s) when computing their difference to lowest order and we end up 
with 



n 



9) 



{2hi 



n 



'l + g)w{g), (37) 



where Ni = q + qf — gqi^ N2 = q + q2 — gqf- Although the integrand yielding 
w{g) is antisymmetric under qi ^ q2, the integral is divergent by power counting 
for n = 2, i.e. it must be computed in = 2 — e dimensions. Its derivative, w'{g), is, 
however, convergent by power counting and it has also an antisymmetric integrand, 
therefore this latter may be evaluated in n = 2 dimensions giving w'{g) = 0. Then 
to compute w{g) one may choose e.g. the point g = —1: 



d^'q qf - ql 
(27r)" (g2 + 1)2 



n 



d'^q 



q^ 



{2nY (g2 + 1)2 



where writing the second equality, we used that gi is n — 1 dimensional, while q2 
is a 1 dimensional variable. From (|3^) one finds that after taking into account the 
change of the renormalization scheme ([23|), in PT the free energy densities of the two 
models (|I]) and (^) do indeed coincide for the LEFT case. Recently this calculation 



has been performed also at the two-loop level fl^ . 
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To discuss the RIGHT and DIAG cases we find it more convenient to use the 
parametrization of for the SU(2) valued field, G\ 

where \1/ resp. $ is a complex resp. a real scalar field. Now ?7r(1) acts as a shift, 
$ I— s> $ + K, and then the corresponding Noether charge of the deformed cr-model is 
simply Qr = / dxp^. Qi^ is slightly more complicated when expressed in terms of 
the canonical momenta (as e^'^" ia"^ = ■ia^e"*'^'^ ): Ql = / c?x[p$ + 2i(p^\l/ — . 



Using Buscher's rule |jTO|, the Lagrangian of the axion (dual) model now takes 
the form: 



^ -soTTV '^^ ^x^^^^ ' ~2' (f^^-f^h (40) 

where / is the dual to $, and 



i^ = ^^-9^^*-^*9^^ = ^(^5^^*-MM^), iV=l + |^|'. (41) 
The canonical transformation connecting Cj] and C^, maps and Ql to 

Qr = -J dxf, Ql = J dx[-f + 2t{p^^-p^*^*)] , (42) 

i.e. Qr Ql do indeed contain the topological charge of the axion field (propor- 
tional to /). 



Applying now the general framework, Eqs. (^^-|33|) to the present cases; i =(L,R), 
= (/,^,^*), with 6*" = (/), y'' = (^, ^*). Using the explicit form of pq, and pq,* 
one finds for J^' ^: 



C.^-O, B-^\ (43) 



C\ - ( Ma B^^J-I^' -fjfl: ^ - f (44) 

where Af = - 2(1 + ^)|^|2)/iV. Substituting these C\ and B'^ into Eq. (H) 

reveals that Hes depends only on and that its extremum is at \l/ = = In 
the DIAG case the actual value of the ground state energy density at this extremum 
is given by: 

iy|min = -^2^(/^R + /^L)', (45) 

A 
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and the expectation value of /' is: /q = —4(1 + g){h-^ + hi^)/X. We note that -ff|min 
agrees (as it should) with the corresponding result of the deformed a-model with 
X i-^ X, g i-^ g (Eq. (3.20) in p). For the RIGHT case the analogous expressions of 



the axion model are simply obtained from (45) by setting /il = 



At this point we recall the somewhat unusual feature of the axion model once 
more, i.e. that the U(l) symmetry of the S'-matrix (^ is realized through a topologi- 
cal current analogously to the Sine-Gordon theory. To emphasize this we quote here 
the value of the classical free energy density corresponding to the Noether charge 
of the 'manifest' Ue(l) symmetry of the Lagrangian (|^): -ff||nin = ~2/i|^/A, quite 
different from Eq. (^5)) with /il = 0. 



To obtain the one loop correction to the free energy, one has to expand C'^ 
around the minimum, Eq. (|45|). Writing / = x/q + /, where the expectation value 
of / vanishes, one finds that the quadratic terms containing / are independent of 
/iL, /iR, so the quadratic pieces of (hence the one loop correction) are effectively 



determined by ^ = y X/2ijj only: 

£2 = d^-^df^ifj* - + «g(^V'* - ^*^) , (46) 

where = 4/iL(/iL^o+^R(l+^o)) and q = /il(1— ^o)~^r(1+^o)- After continuation 
to Euclidean space and writing ip = (0i -t- i02)/\/2, Eq. (|46|) becomes identical to 
the corresponding pieces for the deformed cx-model, Eq. (C.ll) in p (with Aq ^ Xq, 
go ^ go)- From this it follows that the free energy densities fully agree also for the 
DIAG and RIGHT cases in both models. 

Now the results of the comparison between the free energies computed by the 
TEA (based on the proposed S matrix (^) and in PT for the deformed a- model 
@ in Ref. ^ can be simply taken over for the axion model. The conclusion is that 
there is complete consistency between the TEA and the perturbative calculations, 
providing good evidence for the validity of the proposed S matrix (^) for the axion 
model. Since the effective coupling ( pSj ) is identical in the two models the m/Ajj^ 
ratio found in stays unchanged. 

Finally we would like to point out that it would be interesting to study the axion 
model by lattice Monte-Carlo simulations. This would provide us with a completely 
non-perturbative way of testing quantum T-duality. Technical difficulties arising 
from the non-reality of the Euclidean action in the context of the lattice Monte- 
Carlo study of the 0(3) model with a constant 6 term and a suggestion how to 
circumvent them is discussed in |1^ . 
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